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Gravitational Waves from Numerical Mergers 


NR black-hole merger simulations produce waveforms decomposed into (spin-weighted) 
spherical harmonics: 

rMt , r, 9,<t>) = Y^ C em (t, r) _ 2 Y e m (0, <j>). 

£m 


9 We work with strain-rate h = f ipf dt. 

9 Each mode has an amplitude and complex 
phase: 

rh tm = A tm (t)er*«-V) m 

9 A handful of modes dominate energy flux; 
mostly (l, ±1). 

9 (2, ±2) is sufficient for detection ; other 
modes are important for parameter 
estimation. 
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Implicit Rotating Source Picture 


[Baker et al. (2008), Kelly et al. (201 1 )]: 

O Most important WF modes had consistent 
rotational phases m <pem/m 

O Best matches for C = m modes 

9 Rotational frequency model is a smoothed 
“step function” to fundamental QNM 
frequency: 



Fit of fi(f) for (2, 2) mode of xi_uu (top) 
and xi dd (bottom). 
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9 Poor for times earlier than ~ 60M before 
merger 

O Still need model for individual mode 
amplitudes 
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Aligned Spins: Frequency Evolution 


a Asymptotes to quasinormal mode 
frequency for (£, ±t) 

a X4_00 tracks with xi_00 before merger, 
xi_dd after 

a How do x4_oo & xi_dd differ? 

O All BHB mergers result in Kerr hole 
parametrised by mass /W f , spin 

3f = Mfjf = Sf/ M{. 

a Hints at simplest way of breaking QNM 
degeneracy in waveform models 

a Using end-state predictions for 
aligned-spin mergers, can construct 
“equivalence classes” of BHBs 



-60 -40 -20 0 20 40 

t ~ £peak(2,2)(MAij) 

uji m for m = £ modes. 



Equivalence Classes of BHB Mergers 


Aligned-spin BHB configurations “equivalent” to nonspinning 1 :1 , 4:1 , and 6:1 BHBs. 


run name 

q = /M 2 

h 

k 

k 

X1_00 

1.0 

0.0 

0.0 

0.6865 

X1_UD 

1.0 

0.8 

-0.8 

0.6865 

X5_U0 

5.0 

0.0937 

0.0 

0.4748 

O 

o 

1 

X 

4.0 

0.0 

0.0 

0.4748 

X3_D0 

3.0 

-0.1557 

0.0 

0.4760 

X2_D0 

2.0 

-0.4734 

0.0 

0.4765 

X1_DD 

1.0 

-0.6365 

-0.6365 

0.4765 

o 

o 

1 

X 

6.0 

0.0 

0.0 

0.3762 

X5_D0 

5.0 

-0.0752 

0.0 

0.3762 

X4_D0 

4.0 

-0.1896 

0.0 

0.3762 

X3_D0 

3.0 

-0.3843 

0.0 

0.3762 

X2_DD 

2.0 

-0.6250 

-0.6250 

0.3762 
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Modeling IRS Multimode Amplitudes 


IRS model [Kelly et al. (201 1 )] suggested a 
general form for IRS amplitude functions: 


^m(f) — ^0 - 


Ht) 


\| 1 +cm { f 2 - f A )' 


Three free parameters for each (£, m) pair: A 0 , 
n, a -\ . 





for 4:1 nonspinning merger. 
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Modeling IRS Multimode Amplitude Parameters 


® Assemble broad set of aligned-spin BHB 
mergers. 



h 



ii 


BHB configurations used. 
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Modeling IRS Multimode Amplitude Parameters 


9 Assemble broad set of aligned-spin BHB 
mergers. 

O Collect amplitude params {A 0 , n, a } over 
BHB configurations. 



A 0 for (2,2) mode of all configurations. 
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Modeling IRS Multimode Amplitude Parameters 


9 Assemble broad set of aligned-spin BHB 
mergers. 

O Collect amplitude params {A 0 , n, «•, } over 
BHB configurations. 

O Easier to model cuts along BH parameter 
directions . . . 

O Symmetric mass ratio 

ri = Mi M 2 /(M 1 + M 2 f < 0.25 


A 0 for (2,2) mode of nonspinning 
configurations. 
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Modeling IRS Multimode Amplitude Parameters 


9 Assemble broad set of aligned-spin BHB 
mergers. 

O Collect amplitude params (A 0 , n, «•, } over 
BHB configurations. 

O Easier to model cuts along BH parameter 
directions . . . 

O Symmetric mass ratio 

r) = Mi M 2 /(M 1 + M 2 f < 0.25 

O “Total” spin / = (q 2 ^ + k)/(q 2 + 1) 



and 3:1 mergers. 
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Modeling IRS Multimode Amplitude Parameters 


9 Assemble broad set of aligned-spin BHB 
mergers. 

O Collect amplitude params {A 0 , n, «•, } over 
BHB configurations. 

O Easier to model cuts along BH parameter 
directions . . . 

O Symmetric mass ratio 

ijsMt M 2 /(M 1 + M z f < 0.25 

O “Total” spin / = (q 2 ^ + j 2 )/{q 2 + 1) 

O Simplest fit model is product of mass-ratio 
and spin forms: 

MvJ) = g{v)-h(j) 
g(v) = + + 

h(j) = 1 + hi] + h 2 ~j 2 . 



Fit to A 0 for (2,2) mode of all mergers. 
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Modeling the (3,2) Mode 

Kelly & Baker (2013) showed that “observed” (3, 2) mode at merger is largely (2, 2) mode, leaked 
through mismatch between spherical harmonics _2 and spheroidal harmonics _ 2 yf. 


O Visible in both amplitude and frequency. 



-60 -40 -20 0 20 40 60 


t ~ tpeak(2,2)(^AH) 

Amplitude and frequency for several runs. 
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Modeling the (3,2) Mode 


Kelly & Baker (2013) showed that “observed” (3, 2) mode at merger is largely (2, 2) mode, leaked 
through mismatch between spherical harmonics _2 and spheroidal harmonics _ 2 yf. 


9 Visible in both amplitude and frequency. 

9 QNM eigenfunctions are spheroidal 
harmonics 

O Overlap with spherical harmonics is 

s l'lm = dQ _ 2 ^ m (a f (j 22 ; 9, <p) _ 2 Y^{9, <f>)* 


9 .. . leading to mixing coefficients 


Pbasis,^2 — 


Sl'22 

S 2 / 2 2 



Pbasis,» over range of final Kerr spins. 



Modeling the (3,2) Mode 

Kelly & Baker (2013) showed that “observed” (3, 2) mode at merger is largely (2, 2) mode, leaked 
through mismatch between spherical harmonics and spheroidal harmonics _ 2 iy™. 


a Visible in both amplitude and frequency. 


a QNM eigenfunctions are spheroidal 
harmonics 

O Overlap with spherical harmonics is 

seem = j c/n _ 2 ^ m (a f <7 22 ; 6 , <j>) _ 2 <f>Y 


a . . . leading to mixing coefficients 


Pbasis,£2 — 


Sl'22 

S 2 '22 


a Fits observed (3, 2) modes very well; 
(4, 2) numerics too uncertain 



□ ► < g ► <!► <!► 1 '0 ( \0 


APS 2013 8/11 


Kelly (NASA GSFC) 


IRS Multi-Mode Waveforms 



Debumping the (3,2) Mode 


In principle, we can now take the observed (3, 2) mode (red), subtract (2, 2) leakage (blue) to 
yield beautifully single-peaked “real” (3, 2) modes (green). 
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Debumping the (3,2) Mode 


In principle, we can now take the observed (3, 2) mode (red), subtract (2, 2) leakage (blue) to 
yield beautifully single-peaked “real” (3, 2) modes (green). 




Debumping the (3,2) Mode 


In principle, we can now take the observed (3, 2) mode (red), subtract (2, 2) leakage (blue) to 
yield beautifully single-peaked “real” (3, 2) modes (green). 



. . . not as beautiful as hoped. Other mixing effects may still be important before merger. 



Summary 


9 IRS picture gives useful phase model for late-merger and hybrid waveforms. 

9 Individual mode amplitudes require separate parameter models (as in post-Newtonian). 
9 (£,m< £) modes need to have “leakage” subtracted before fitting. 

9 But more subtraction might be necessary . . . 

9 Next: (1) test against ringdown templates for detection & parameter estimation. 

9 Next: (2) use higher-accuracy WFs from the NR-AR collaboration to generate better fits. 
9 (or CCC catalog . . . ) 
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